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Entropy and the fluctuation-dissipation theorem are at the heart of statistical mechan- 
ics near equilibrium. Driving a system beyond the linear response regime leads to (i) the 
breakdown of the fluctuation-dissipation theorem and (ii) a nonzero entropy production 
rate. We show how both phenomena are related using the general framework of stochas- 
tic thermodynamics suitable for soft matter systems governed by stochastic dynamics and 
driven through nonconservative forces or external flows. In particular, the excess of the 
fluctuation-dissipation theorem in a nonequilibrium steady state compared to equilibrium is 
related to total entropy production. Alternative recent derivations of generalized fluctuation- 
dissipation theorems are sketched and related to each other. The theory is illustrated for 
two systems: a driven single colloidal particle and systems driven through simple shear flow. 



§1. Introduction 

The fluctuation-dissipation theorem (FDT) is one of the cornerstones of statis- 
tical mechanics. Going beyond equilibrium, it connects the response of a system 
perturbed slightly out of equilibrium with correlations in equilibrium lif Specifically 
for a system at temperature T, and Boltzmann's constant set to unity, the FDT 
reads 

-d t (Mt)B^\t')) , (i-i) 



TR ^ {t ~ t)= ^h(¥Y 



h=0 



i.e., the response of the system after some perturbation h has been applied is mea- 
sured through a time-dependent change of the mean value of A. This response is 
equal to the time-derivative of a correlation function in equilibrium involving A and 
another observable B^ eq \ This second observable at earlier time is not arbitrary but 
rather is the conjugate of the field h with respect to the system's energy U, i.e., ap- 
plying the field changes the energy as U ^ U-hB^. Onsa ger's regression principle 
casts this remarkable symmetry into words: the decay of spontaneous fluctuations 
cannot be distinguished from the decay of a forced fluctuation. 



The values of h for which the form (IT) of the FDT is valid determine the linear 
response regime. Concepts like linear irreversible thermodynamics^! 1 and local equi- 
librium assumptions have been used successfully to extend certain thermodynamic 
concepts into this near-equilibrium regime based on the notion of entropy produc- 
tion. However, a general theory for going beyond the linear response regime is still 
missing. In practical terms, the absence of a concept comparable in universality 
to the Gibbs-Boltzmann distribution is probably the biggest obstacle in formulat- 
ing and applying a general nonequilibrium thermodynamics. During the last 10-15 
years substantial progress has been made on another front through the formulation 
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and study of nonequilibrium fluctuation relations valid arbitrarily far from equilib- 
rium. These relations constrain the probability distributions of quantities like work, 
heat, and entropy production. The arguably most famous representatives are the 
nonequilibrium work relations due to Jarz' and CrookaHriHP and the fluctu- 

ation theorems for entropy productionEfli^l Subsequently, it has been shown that 
thermodynamics can be formulated consistently for driven systems on the level of 
single trajectoriesE3~fHl' 

The systems of interest to this work are soft matter systems such as colloidal 
suspensions, single colloidal particles, biomolecules such as DNA and RNA, as well 
as motor proteins such as Fi-ATPase. These systems share the property that they 
are immersed into a host fluid of well defined temperature and that they intrinsically 
operate in nonequilibrium; they are either driven by chemical gradients, mechanical 
forces, or external flows. However, we will require the fluctuations arising from the 
bath-system interactions to be described by equilibrium fluctuations. This cannot 
be strictly true, for an explicit calculation see, e.g., RefJH However, several exper- 
imental tests for single colloidal particleaiSlllZf confirm the theoretical predictions 
for measured probability distributions, thus supporting the validity of this approxi- 
mation. 

The FDT is one of the ubiquitous tools in statistical mechanics and computa- 
tional physics. Due to its importance, possible extensions into the realm of nonequi- 
librium, especially for glassy dynamics, of been studied for a long time leading to 
a number of reviews (a selection is Refs.jHf^HH) see also references therein). In this 
paper, based on Ref.pU we mainly discuss the connection between entropy produc- 
tion as defined in the framework of stochastic thermodynamics and the fluctuation- 
dissipation theorem. The main ingredient is that, although the system is driven 
beyond the linear response regime, a linear response in reaction to a slight perturba- 
tion of the driven system out of its nonequilibrium steady state can still be defined. 
Hence, while the equilibrium fluctuation-dissipation theorem no longer holds in the 
form (1-1), it can be extended to nonequilibrium steady states by an additive cor- 
rection. The principle of conjugate observables with respect to energy is replaced by 
conjugate observables with respect to entropy production. We discuss the relation 
to other recent work on the FDT and finally we give specific expressions for the case 
of a single colloidal particle and colloidal suspensions or polymers driven by simple 
shear flow. 



§2. Driven Soft Matter 



In the following, we consider stochastic systems obeying Markovian dynamics. 
The state space might either be discrete or continuous, and x denotes a single ele- 
ment in this space. For time-independent forces or rates, the system will eventually 
settle in a steady state with probability distribution ipo(x; {Xj}). It is important 
to note that this steady state depends on parameters Xj that we control externally. 
To simplify notations, however, we will often not write this dependence explicitly. 
The values of these parameters determine whether the steady state is equilibrium 
or a nonequilibrium steady state characterized by a non-vanishing mean entropy 
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production rate. 

2.1. Continuous state space 

Soft matter systems such as colloidal particles or biomolecules immersed in a 
fluid are often well described by overdamped Langevin dynamics. The fluid acts 
as the heat reservoir, i.e., we assume that the fluid remains in equilibrium with 
temperature T even though we drive the immersed subsystem. In addition, the fluid 
might have an imposed flow profile u(r) . We consider the system to be composed of iV 
'units' (colloidal particles or monomers) with positions rjt, ignoring internal degrees 
of freedom. The configuration of the system is then given by x = {r±, . . . , r^r}. The 
time-evolution of the probability distribution ip(x,t) is governed by the Smoluchowski 
equation 

N 

d t iP + Y,V k -{vki>) = 0, (2-1) 

k=l 

where 

v fe = u(r fc ) + fi F k = u(r k ) + /i [-V k U + f k - TV k In %b] (2-2) 

is the local mean velocity and fiQ is the bare mobility due to friction. Any deviation 
of the local mean velocity v*. from the imposed flow profile u(r) has to be caused by 
a force F k exerted on the fc-th particle. We allow for three contributions to this force: 
(i) conservative forces S7 k U due to the potential energy U(x), (ii) nonconservative 
forces ffc, and (iii) "thermodynamic" forces arising from the stochastic interactions 
between system and the surrounding fluid. In equilibrium, i.e., in the absence of 
external flows and nonconservative forces, detailed balance holds which amounts to 
Vfc = 0. The thermodynamic forces ensure th at th e equilibrium Gibbs-Boltzmann 



distribution ipQ oc e~ u ^ T is the solution of Eq. (2-1). Although a force on a particle 
in principle leads to a distortion of the flow and the coupling of forces, here we will 
ignore hydrodynamic interactions between particles. 

An equivalent dynamic prescription on the level of single stochastic trajectories 
is the Langevin equation 

h = to [-VfcU + f fc ] + £ fc , (2-3) 

where the noise £ k explicitly models the stochastic interactions between particles 
and the surrounding fluid. The noise has zero mean and correlations = 
2fi Tt5 kk/ 5{t-t'). 

2.2. Stochastic thermodynamics 

Stochastic thermodynamics is a conceptual framework combining energetics along 
single trajectories with the definition of a stochastic entropy We first extend the 
idea by SekimotdHJ'lllr to define work and heat along single stochastic trajectories 
to the situation where external flows are present Hlf In a thermodynamic context, 
the work is the change in energy that is controlled externally. The work rate readelll 

dU N N 



J 



k=l k=l 
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u=-v 



Fig. 1. Colloidal particle dragged by optical tweezers through a viscous fluid, a) Laboratory frame 
of reference, b) Co-moving frame of reference: an observer moving such that the trap appears 
stationary would see a flow trying to advect the particle. The work is invariant in both scenarios. 

The first term is the work spent to change the potential energy through changing 
control parameters. The other two terms are due to the external flow and the 
nonconservative forces, respectively. Conservation of energy in form of the first law 
of thermodynamics then leads to the heat rate 

dU N 

q = w - — = J2i~ V kU + f fe ] • [r k - u(r fe )]. (2-5) 
k=i 

Hence, the total force times displacement of the particles equals the heat dissipated 
into the surrounding fluid. 

To make these ideas transparent and as an example, consider the prototypical 
colloidal particle dragged by optical tweezers, see Fig. [T] The optical trap is modeled 
as a harmonic potential with strength k and focus position A. This system has 
been studied comprehensively both theoreticall^HH 1 'IHI and experimentally The 
parameter we control externally is the position of the trap center. Alternatively, one 
could control the trap strength. The potential energy reads U(x; A) = (k/2)(x — A) 2 



and the work rate following Eq. (2-4) is w = Xk(x — A), where x is the position of the 
particle in the laboratory frame. An observer moving with (constant) velocity v = A 
measures the particle position y = x — vt. As expected, the expression for the work 
rate, w = —kyv, is invariant under such a change of the frame of reference although 
A = since now the fluid appears to flow with velocity u = —v. The corresponding 
expressions for the heat are q = —k(x — X)x = —k{v + y). 

The final step is the introduction of a stochastic entropy defined a 

S (t) = -ln^(x(t),t;{A,(t)}). (2-6) 

Averaging this expression leads to the well known Shannon entropy 

(s(t)}=- dxip(x,t)lnip(x,t). (2-7) 



Taking the time-derivative of Eq. (2-6), we obtain the balance equation 

ds 

= — Smed + «tot, (2-f 
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i.e., the total entropy production rate s t ot is the sum of the change of entropy of the 
system ds/dt and the change of medium entropy s me d- Since we demand the fluid 
to stay in equilibrium, we can identify the dissipated heat with the change of fluid 
entropy through the Clausius relation s me d = q/T. 

2.3. Discrete state space 

Discrete states might be intrinsic but often they arise through some spatial or 
temporal coarse graining procedure. A biochemically motivated example for the for- 
mer is the discrete motion of Fi-ATPaseHU Transitions between states occur with 
rates W(x — > y). Analogous to the continuous case discussed in the previous section, 
the thermodynamic notions of work and heat can be defined consistently H^f How- 
ever, stochastic entropy itself might be introduced in a more abstract way regardless 
of whether the system of interest is coupled to a heat bath© A single stochastic 
trajectory x(t) of length i Q b s consists of K jumps at times < t a < t Q ^ s from state 
x a -i to state x a with a given initial state xq. The total change of stochastic entropy 
along this trajectory (in a steady state) can be written 

As ^ s(t ohB ) - s(Q) = MXa) 



K K 

E W(x a -> X a -i) _ 1po(x a )W(x a ->■ X a -1 



^ W(x a -x -)• x a ) ^ i/j (x a -i)W(x a -i -> x a ) 

= -S me d[x{t)] + Stot [»(*)], 

(2-9) 

where in the second line we have expanded the sum by inserting the transition rates. 



The sums run over all transitions. Following Eq. (2-8), we interpret the first sum as 
the change of medium entropy with rate 

s mcd (t) = t a ) ln g^LZ^g=0 (2-10) 

, W(X a -i ->■ X a ) 



The second term in Eq. (2-9) is the total entropy production. In equilibrium, detailed 
balance holds with i/jq(x)W(x — > y) = ipo(y)W(y — > x). Hence, in equilibrium the 
total entropy production becomes identically zero. Moreover, (stot) ^ 0, where the 
equal sign holds in equilibrium. This identification of entropy and the corresponding 
fluctuation relations for their probability distributions have been illustrated experi- 
mentally for a driven single defect in diamond j22f'H3 



§3. The Fluctuation-Dissipation Theorem 



Before discussing the FDT in nonequilibrium, we will briefly introduce an ab- 
stract notation particularly apt to treat continuous and discrete state space on equal 
footing and to see the general structure of the FDT we want to unveil. Using the 
bra-ket notation, the system is described by the state vector \ip(t)) from which the 
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probability distribution is read off as ip(x,t) = (x\ip(t)). We define a reference state 
(— | through (—\ip(t)) = 1 which expresses the conservation of probability. The time 
evolution of the state vector is given through 

dtm)) = L\m) (3-i) 

with an operator L obeying (— \L = 0. The steady state \ipo) corresponds to the 
right eigenvector, L\^o) = 0, of this time evolution operators with eigenvalue zero. 
In this picture, the observables correspond to operators with expectation values 

(A(t)) (-\A\m)- (3-2) 

For diagonal observables, (x|A|a/) = A(x)5(x — x'). This expression reduces to the 
familiar {A(t)) = J dx A(x)^(x,t) upon inserting the identity 1 = J dx \x){x\. 
Similar expressions where the integration is replaced by a sum are obtained for 
discrete state spaces. We denote the average value of A in the steady state by (^4)o- 
Finally, two-point correlations of two observables A and B in the steady state become 

(A(t)B(0)) o = (-\Ae tl B\^ ). (3-3) 

For completeness, we give two explicit expressions for L. First, the matrix elements 
of the time-evolution operator for the Smoluchowski equation (2-1) read 

N 

(x\L\x') = - Vfc ■ {u(r fc ) + Mo[-V fc *7 + f k - TV k }} 5(x - x'). (3-4) 

k=l 

Second, for a general Markov process on a discrete state space, the time-evolution 
operator is the usual left stochastic matrix with elements 

(x\L\x') = W(x' ^x)- 5 XX > W(x -»• y). (3-5) 

3.1. Formal derivation 

Let A be one control parameter, and without loss of generality let A = for the 
steady state. Changing A will perturb the system away from its steady state with 
the following evolution governed by Eq. (3-1). For small A, the system will response 
linearly in the sense 

(A(t)) = (A) + f dt' R AX (t-t')X(t') + 0(X 2 ), (3-6) 

J — oo 

i.e., the mean at later times is a linear functional of the perturbation A(i). The 
response function RA,\(t) connecting both is formally given through 



d\{t> 

Due to causality, the response function vanishes for t! > t 



(S < *)• (3-7) 

A=0 
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The calculation of the response function (3-7) for general Markovian processes 



is well known. For completeness and later reference, we repeat the derivation here, 
followi ng th e route of AgarwaMI and Hanggi and ThomasHH The formal solution 
of Eq. (3-1) for a, due to the perturbation, time-dependent L(t) reads 



W)) 



f-^t' L(t>) 



1*00 



(3-8) 



given that the system at t — > — oo has been prepared in the steady state. The 
exponential is to be understood in the time-ordered sense. W e the n split the time 
evolution operator, L(t) = L + X(t)SL, and use the expression (3-2) to obtain 



= (-\Ae^ L 5L\^) = (A(t)B(t')) 



(3-9) 



after performing the functional derivative with respect to A(i'). Hence, the response 
function can be expressed as a correlation function of the observable A with the 
operator 5L. Although formally correct, it is of little practical use as long as we do 
not find a 'physical' representation B of this observable in the sense that it can be 
expressed in terms of, in principle, measurable quantities. For a first representation, 
we assume the observable 5L to be diagonal, 5L = J dx \x)B^(x)(x\, with 



B^\x) 



{x\SL\j; } 
(x\tpo) 



(3-10) 



This result indeed has been known for more than thirty years and we call it the 
'Agarwal' form. However, B^ still seems to have no transparent physical mean- 
ing. Moreover, for complex systems with a large number of degrees of freedom, 
the explicit form of the stationary distribution ipo{x) in general is neither available 
experimentally nor from numerical simulations. 

In the next subsections we will discuss two more representations of B labeled by 
different superscripts J3W. To this end it is crucial to realize that a whole class of 



representations for the observable B exist that all lead to the same FDT (3-9). The 



formal reason is that we have some freedom in expressing the state 5L\ipQ) . Beyond 
the representations discussed below, there are in principle infinitely many variants 
of the FDT since with B^ = B^ 2 \ where = denotes the equivalence of observables, 
any normalized linear combination 



1 



ci + c 2 

with ci t 2 real will be admissible. 
3.2. Role of stochastic entropy 



c x B {1) + c 2 B {2) 



(3-11) 



We want to establish the connection between the stochastic entropy (2-6) and 
the FDTH3 We consider two steady states separated by a small A with state vectors 
\ipo) and |V>o) + A|(fa/>), respectively. Using that these state vectors are the right 
eigenvectors of the corresponding evolution operators with eigenvalue zero, 



L\6i>) = -SL\^ 



(3-12) 
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holds to first order in A. Inserting this into Eq. (3-9), we obtain 







R A ,x(t - = ^(-\Ae^ L \5^ = (A(t)B^(t')) 
which is independent of 8L. The physical observable is now given as 



£(e), 



X, X 



d dtp(x) 
dt ip (x) 



d ds 
dXdi 



(3-13) 



(3-14) 



A=0 



The total time derivative in the first expression is along a single trajectory starting 
at x. One should keep in mind that this expression is to be averaged over trajectories 



in the correlation function (3T6). Expanding Eq. (2-6) up to first order in A leads to 
—5il>(x)/'ipo(x) a nd, fi nally, we interchange the order of derivations to obtain B^ e \ In 
contrast to Eq. (3-10), B^ in general depends also on x and hence is a non-diagonal 
observable. 



The conceptual advantage of Eq. (3-14) is that it leads to an observable which 



is conjugate to the perturbation parameter with respect to entropy production, in 
the same spirit that the observable £>( cq ) is conjugate with respect to energy in 
equilibrium. To recover the equilibrium form (1-1) of the FDT, we employ the 



stationary distribution in equilibrium given by the Gibbs-Boltzmann distribution 



-[U(x;\)-F(\)]/T 



(3-15) 



where J- is the free energy. The system entropy simply becomes Ts(t) = U (x(t)) —T. 



Eq. (3T3) then reads 



TR^(t-t') 



^(A(t)[-d x U] x=0 (t')) + -^ 



{Ay 



AT 
'dA 



(3-16) 



The free energy is constant along trajectories. The term in the square brackets 
is then constant and the second term vanishes. The expression d\s = (l/T)d\U 
appears in the FDT, which acquires the well-known form (IT). Hence, 
identify B^(x) -- 



we can 



d\U(x; A)| A=0 as expected. 



Considering the general case of a nonequilibrium steady state, the observable ( 3T4 ) 



can be used to find a link to the equilibrium case (IT). Using Eq. (2-8), we split 



5 (e) = d\s med - d\s t 



ot 



(3T7) 



into two terms. In it has been shown that the second term vanishes in 

equilibrium. Hence, in equilibrium the FDT (3-9) can then also be written in the 
form 



R { 2 q ht-t') = (A(t)[d x s med ](t')} . 



(3T8) 



Now consider the same system driven into a nonequilibrium steady state by a driving 
force corresponding to the parameter A. Perturbing the system by a small change of 
the same driving force will leave d\s me d unaltered. We can thus keep the correlation 



function on the right hand side of Eq. (3T8), now evaluated under nonequilibrium 
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conditions, and subtract the second term in Eq. (3-17) that involves the observable 
conjugate to total entropy production to obtain 

RA,x(t - = (A(t)[d x s mcd ](t')} - (A(t)[d x s t ot](t'))o. (3-19) 

The second correlation function quantifies the excess compared to equilibrium. Such 
a splitting has been mentioned befordHMHJ without recognizing the meaning of the 
excess term as a correlation with the total entropy production. Harada and Sasa 
have discussed the connection of this excess with energy dissipation!!!! 

3.3. Path weight approach 



An alternative way to write the average (3-2 ) is to use the path integral formalism 
leading to 

r fx(t)=x 
(A{t)) = / dx A(x) / [dx(t)] P[x(t); X(t)], (3-20) 



where the path integral runs over all trajectories which end in x at time t. The weight 
of a single trajectory is P[x(t)] = e~~ s ^- x ^ with stochastic action S. Calculating the 
response function (3-7) by taking the functional derivative, we obtain 

dxdx' A(x) [dx(t)} P [x(t)}B^\x',x')Mx'), (3-21) 

Jx(t')=x> 

where Pq is the path weight in the steady state and the path integral now sums over 
all paths starting in x' at earlier time t' and ending in x at time t. We therefore again 
find the FDT in the general form (|3-9|) but now with the observable representation 

= -d x C(x,±). (3-22) 

A=0 



b <p),„ ^_ SS[x(t);X(t)] 



ix.x) 



5X(t') 



For the last equality, we have used that the action can be written as S = f dt C 
The explicit expressions reacj*^ 

1 N 



with £ k given through the Langevin equation (2-3) and 

A" 



C(t) = W(x(t) -»• y) - J2 S (t ~ ta) ^ ^(x Q _i x a ) (3-24) 



a=l 



for continuous and discrete state space, respectively. 

A similar approach has been used by Baiesi et al. to derive yet another form 
of the FDT by relating the path weight of the perturbed process with the station- 
ary path weightJ!!J However, certain forms of this FDT have been know for a long 



There is a Jakobian involved in the change of variables from £ to x which we dropped since 
it does not contribute to —d\C 
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timell^rl 3 !! For example, realizing that a force perturbation is equivalent to a per- 
turbation of the noise leads immediately to 

TR A>fk (t-t') = » (^^ = \{A{t)t k {lf)). (3-25) 

This has been exploited in RefJD 

3.4. Hatano-Sasa relation approach 

The term steady state thermodynamics has been coined for a phenomenological 



theorjc2rl!2f promoting the splitting of the total dissipated heat into a housekeeping 
heat (/hie and an excess heat q ex . We introduce a pseudo-potential 4>{x) via the steady 
state probability 

tfo(x) = e~^ x) , <t>{x) £ U(x)/T (3-26) 

in analogy to, but different from, the Gibbs-Boltzmann distribution. The excess heat 
is q e x/T = Y — Acf), where Y is the transition functional defined as 

y M t )1 = /«A. W *<y>» (3-27) 

For this transition functional, the fluctuation relation 

(e~ Y ) = 1 (3-28) 
holds ED For completeness, note that also the housekeeping heat c/hk = Q~ Qcx 

fulfills 

a fluctuation relation© 
We define the variables 

X, = - * (3,9) 

conjugate to the paramet er A.,- with respect to the pseudo-potential <f>. Expanding 
the Hatano-Sasa relation (3-28) in powers of the A's, Prost et a/. 431 derived the FDT 

(Xj) = I dt' R jk (t - t')X k (t'), R jk (t - t') = -5 t (X i (t)X fe (i'))o, (3-30) 

J — oo 

where (-Xj)o = by construction. Noting that in a steady state the stochastic 
entropy (2-6) becomes s(t) = <f>(x(t)), we see that this result is equivalent to the 
FDT (3-13). However, Eq. (3-13) seems to be more general since it does not restrict 
the observable that measures the response. 

§4. Illustrations 

For an illustration, we consider two systems that have been discussed in detail 
previously HD'Ell sheared soft matter systems and a single particle moving in a 
periodic potential. 
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4.1. Shear driven systems 

The system consists of N particles with positions x = {ri, . . . , r n |. The FDT 
in such shear driven suspensions has been addressed numerically^ 1 &j and in the 
framework of mode-coupling theory EfMl3 Invariant quantities^} constitute an ex- 
act result for systems driven through boundaries with an otherwise unaltered Hamil- 
tonian. In contrast, in this example the system is driven through an imposed flow 
with profile u(r) = (7j/,0,0) T , where 7 is the strain rate. Shearing the fluid builds 
up stress with off-diagonal element 



N 



0- X y(x) =^Vk 



k=l 



dU 

dx k 



(4-1) 



which will be the crucial quantity in the following. 

The time-evolution operator is given by Eq. (3-4). Perturbing the strain rate 
A h> 7, we obtain 5L = — Ylk=i ^ fc af~- ^he ^ rs ^ representation of the conjugate 
observable is the diagonal 'Agarwal' form (3-10) 

5Lip (x) 



a 



xy\ 



0~xy{%)- 



The second term here is 

N { . 

o- xy (x) = ^2[U(x) - T(j)(x)] = ^Vk{vk 



N 



k=l 



/'0 



IVk) 



(4-2) 



(4-3) 



k=l 



which involves the x-component of the local mean velocity (2-2). The medium en- 
tropy production rate following Eq. (2-5) is 



v 



Ts med = y~][r k - u(r fe )] • {V k U) 



(4-4) 



k=l 



and therefore T9 7 s me d = o~ xy . The representation B^ based on the splitting of the 
entropy production (3-17) is therefore also diagonal and has to coincide with the first 
form, £>( e ) = B^ a \ A different representation is obtained as 



TB (p) 



E- 

r 2 » 



x k + no 



dU 

dx k 



IVk Vk 



1 

Ifxy + 



IVk) (4-5) 



following the path integral approach (3-22) with Eq. (3-23). Combining the two 
different forms, we arrive at 



T[2B^ - B^] 



1 ^ 



Mo 



k=l 



Vk{x k - u k>a 



(4-6) 



This expression also quantifies some kind of stress, where the velocity in x-direction 
is measured with respect to the flow. 
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Fig. 2. Rouse polymer: The integrated response x(i) [Eq. ( 4- 7 1] plotted vs the integrated correlation 
function K(t) [Eq. (4-8 1] for different strain rates 7 and normalized by Koa = K(t — > 00). In 
equilibrium (7 = 0), this plot is a straight line with T\(t) = K(t). 



Knowing the different representations the FDT can acquire might help to mea- 
sure, but of course does not specify, the actual functions R(t), C(t), and I(t). Ex- 
plicit expressions for response and correlation functions have been obtained for 
the case of a Rouse polymeiHU 1 in RefEU 1 The potential energy reads U(x) = 
(k/2)^2 k ~i |r^ + i — rfc| 2 . In the limit N — > 00, the following expressions for in- 
tegrated response and stress auto correlations have been calculated, 



xit) 



dtR a Jt) = T (l-e-*) + \/rterfc(\/i) 



(4-7) 



K(t) = / dt [(a xy (t)a xy (0)) - (a xy ) 2 ] = T 2 (1 + 7 2 )(1 - cT*) + V^ievfc(Vi) . 
Jo L J 

(4-8) 

Here, the time t has been scaled by the fundamental relaxation time T\ = (2fiQk)^ 1 . 
In Fig. [2J the integrated response x{K) as function of the integrated correlation 
is plotted for different strain rates 7. While for 7 = we observe a straight line 
corresponding to Tx(t) = K(t), the excess and therefore the deviation from the 
straight line increase with increasing strain rate 7. 

4.2. Single particle in a periodic potential 

Another example we have studied previously is that of a driven single colloidal 
particle moving in one dimension in a periodic potentialUlMIIf The Langevin equa- 
tion reads 

x = no[-d*U(x) +/]+£ = VoF(x) + e. (4-9) 

Here we choose A 1— > f, i.e., we pe rturb the driving force /. The perturbation 
operator 5L = — fj>od x with Eq. (3T0) produces the diagonal Agarval form B^\x) = 
v{x) — noF(x). Following the path weight approach leads to B^ p \x,x) = (l/2)[i; — 
/xo-F(x)]. A straightforward combination of both results amounts to 



5 (e) = 2£ (p) - £ (a) 



(4-10) 
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Since form Eq. (2-5) TdfS me d = x follows we can identify TdfStot — v, i-e., the 
excess involves as observable the local mean velocity (2-2). The idea of identifying 
the excess with the local mean velocity has lead to a formulation where the FDT 
acquires its equilibrium form through a transformation into the Lagrangian reference 
frame moving with the local mean velocity H°HH Fluctuations in this co-moving 
frame are essentially equilibrium fluctuations and the driving manifests itself in the 
mean values of observables. Of course, this result also crucially depends on fixing 
the bath fluctuations to be equilibrium fluctuations. 



§5. Summary 



In summary, the general FDT for Markovian dynamics is considered. Following 
Ref.jm 1 a generalized FDT is derived where the conjugate observable is determined 
with respect to entropy production of the system in contrast to equilibrium where 
the conjugate observable is determined with respect to energy. Alternative ways to 
derive generalized FDTs for driven systems are sketched based on the path-integral 
and on the Hatano-Sasa relation (3-28). It is shown how the results of different 
approaches can be combined to produce new conjugate observables through linear 
combination. The general theory is illustrated for shear driven systems where the 
strain rate is the perturbation parameter and a driven single colloidal particle with 
the driving force taking the role as the perturbation parameter. 
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